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1 Introduction 



Lattice calculations in QCD have shown significant advances in the last years [1,2]. Sim- 
ulations at the physical value of the pion mass are carried out nowadays, including up, 
down and strange quarks as dynamical degrees of freedom. A further step to render lattice 
QCD calculations even more realistic is also to include the charm quark in the simulations. 
In fact the simulations with two quark generations have already been started [3,4,5,6,7,8]. 
Clearly, having the dynamical charm degree of freedom allows to test many physical as- 
pects of the charm sector of QCD, such as the mesonic and baryonic spectrum and decay 
constants, heavy quark effects in operator matrix elements, the renormalized charm quark 
mass and eventually the running of the strong coupling constant for four flavours. 

One concern when adding a charm quark mass is that lattice spacing effects may become 
large. Even in C(a)-improved lattice theories cut-off effects of 0(a 2 ml h£iTm ) are expected 
to be present which then can become significant due to the rather large value of the charm 
quark mass. In this paper, we want to investigate Wilson twisted mass fermions [9,10,11] 
in a formulation that comprises mass degenerate light up and down quarks and mass 
non-degenerate strange and charm quarks, which we refer to as Nf = 2 + 1 + 1 setup. All 
our calculations were performed at maximal twist, where automatic (9(a)-improvement is 
realized. The particular goal of this paper is to prepare the analytical basis for a study 
of the lattice spacing effects at tree-level of perturbation theory. To this end, following 
Refs. [12,13,14], we derive the quark propagators both for the 4-dimensional discrete 
momentum representation and the time-momentum representation in this Nf = 2 + 1 + 1 
setup. From the expressions of these propagators we construct the full (four times four) 
matrix of the correlation functions for K and D mesons. 

Although calculations at tree-level of perturbation theory cannot be a quantitative mea- 
sure for the interacting case when gluon degrees of freedom are incorporated, the qualita- 
tive results obtained in such setup can nevertheless serve as a valuable indicator. Studying 
the relative size of lattice spacing effects for the pion, the K and D mesons provides impor- 
tant information how cut-off effects behave when light and heavy quarks are considered. 

In the twisted mass formulation of lattice QCD, the situation of maximal twist is of 
particular importance since in this case all physical observables are automatically 0(a)- 
improved. We don't need an operator improvement and the calculation of corresponding 
coefficients. In the interacting case, maximal twist is realized by tuning of the Wilson 
quark mass to its critical value. As numerical and theoretical investigations in the past 
have demonstrated, see Refs. [15,16,17,18,19], there are, however, optimal and non-optimal 
ways to realize maximal twist. 

At tree-level of perturbation theory, we study the analogs of these definitions of maximal 
twist. The optimal definition of maximal twist in the interacting case is to tune the PCAC 
quark mass to zero, which corresponds to setting the bare Wilson quark mass to zero at 
tree-level. The choice of tuning of the theory to maximal twist is not unique. Any definition 
of the critical quark mass that differs from the optimal one by 0(a) effects still leads to 
maximal twist. However, such definitions can introduce unwanted, chirally enhanced 0(a 2 ) 
effects, i.e. terms that go like 0(a 2 /m 2 ) and are hence called non-optimal. Our setup of 
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tree-level of perturbation theory allows us to study both choices, the optimal tuning and 
the non-optimal tuning to maximal twist and we explore both options in this paper to 
understand better the behaviour of meson masses as functions of a 2 , when different tuning 
conditions are employed. 



2 Tree level of QCD with Nf = 2+1 + 1 flavours of quarks 



2.1 Tree level of QCD in the continuum 



2.1.1 Physical basis versus twisted basis 

We consider a light degenerate quark doublet ip^ = {^ u \i/j^} with a mass m u>d and a 
heavy non-degenerate quark doublet ip^ = {ip^ c \ ip^} with masses m c ^ m s . At tree- 
level of perturbation theory, i.e. in the absence of any gauge field, the light and heavy 
quark actions in the physical basis read: 

Si = J d 4 x^ l \x)(j^ + m u>d y^(x), (1) 



S h = J d 4 ^ (ft) (i)( 7 A + diag(m c ,m a ))v (h) (aO- (2) 

Expressing the physical basis quark fields if)^ and ipw in terms of twisted basis quark 
fields an d y i a the twist rotation: 

,0(0 _ e «wnsTs/2^,(0 j ^p{h) _ e iu) h -ysri/2^,(h) ^ 

yields: 

Si = J d 4 x x {l) {x) (j A + m 0>l + i^rzj X (l) (x) , (4) 
with: m ,i = m M cos(w ; ), pi q = m M sin(c^) and 

S h = J d 4 x x (h) (x) f 7 M fy + m 0<h + ifi^n + fi 5 r 3 ^ X (h) (x), (5) 

where: m 0th = ((m c + m s )/2) cos(u h ), \i a = ((m c + m s )/2) sm(u h ) and /i s = (m c - m s )/2 
respectively (the Dirac matrices 7 M [pi = 0, 1, 2, 3] and 75 act in spinor space and the Pauli 
matrices r a [a = 1,2,3] act in the flavour space), ui and Uh are the so-called light and 
heavy twist angles. 

The above relations between the quark masses and the mass parameters in the twisted 
basis quark actions can be solved with respect to the quark masses: 

m u>d = \j{m ,i) 2 + 9 ) 2 , (6) 
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m s = yj (m 0ifc ) 2 + (/i CT ) 2 - fig, m c = yj (m 0>h ) 2 + (fx a ) 2 + fj, s . (7) 

At maximal twist ui = Uh = vr/2, corresponding to m 0i ; = m 0i h = 0, these relations 
simplify to: 

m u>d = pLq, m s = jJL a - jig, m c = fi a + /j, s . (8) 



2.1.2 Meson spectrum 

At tree-level of QCD, i.e. in absence of gluonic fields and, therefore, of any interactions 
between quarks, mesons correspond to free quark- ant iquark pairs. For mesons at rest, i.e. 
with total momentum P = 0, the quark and the antiquark have to have opposite momenta 
+p and —p. For example the spectra of K and D mesons are given by: 



m K (p) = \/ (rn s ) 2 + p 2 + \J (m u4 ) 2 + p 2 , m D (p) = \J (m c ) 2 + p 2 + \J (m Ujd ) 2 + p 2 . (9) 

When we consider infinitely extended space, there is a continuum of states parametrized 
by momentum ±p associated with the quark-antiquark pair. When considering a finite 
spatial volume L 3 , only discrete momenta p = 27rn/L, neZ 3 are possible, rendering the 
meson spectra also discrete. In both cases all states are two-fold degenerate, corresponding 
to parity V = — and V = +. An exception are mesons with quark momentum p = 0, for 
which one can show that at tree-level they only exist for negative parity. 

2.2 Tree level of perturbation theory of twisted mass lattice QCD 

The lattice discretizations of the twisted basis quark actions (4) and (5) are: [9,10] 

Si = a^x^W^wKd + WaJx^W, (10) 

Sh = a 4 Y<X {h) (x)(D w (m 0ih ) + wy 5 7i + W r 3 ))( W W, (11) 
where Dyj denotes the standard Wilson Dirac operator: 

^wK) = \ (V„ + V*) - aV^V M ) + mo,*, x E {I, h}. (12) 



At maximal twist, physical observables are automatically 0(a) improved, i.e. lattice dis- 
cretization effects appear only quadratically [10,11]. Although maximal twist can be real- 
ized in many different ways, there is an optimal definition of maximal twist corresponding 
to setting m 0) ; = mo^ = [16]. Note that for the light quark doublet, the Wilson term 
explicitly breaks isospin and parity, which becomes clear after a rotation to the physical 
basis. Only parity combined with light flavour exchange remains a symmetry. Since the 
Wilson term is 0(a), isospin and parity are restored, when approaching the continuum 
limit. For the heavy quark doublet similar statements apply. 
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For infinite temporal extension, the lattice meson spectrum is expected to be qualitatively 
identical to the continuum meson spectrum discussed in Section 2.1.2. In particular, also 
on the lattice no positive parity mesons exist, with both quarks having zero momentum. 



3 Twisted mass lattice quark propagators 



The calculation of the heavy strange and charm quark propagators at tree- level of twisted 
mass lattice QCD is somewhat lengthy, but straightforward. Here we only quote the result, 
while details regarding the calculation can be found in App. A. 

We can express the twisted mass action of the heavy doublet (11) in terms of the matrix 
K(x;y): 



S h = a*YlY,X ih \x)K(x;y)x W (y), 



(13) 



x y 

which is of the form: 



K(x; y) = J] (Sx+fry (1 - 7 M ) + 5 x _^ y (l + 7^ + 5 XjV ^ (m 0}h + ^ + i/vysTL + figr^j , 



(14) 



where x, y G Z denote space-time indices. The equation: 



(15) 



relates K(x; y) to the heavy twisted mass propag ator SW(y; z) in position space repre- 
sentation. 

In the time-momentum space representation, defined by: 



SW(*,p) = £ e" ip ( x ^SW(x ,x;y ,y), (16) 
where t = Xq — yo, the resulting propagator for infinite temporal lattice extension reads: 



SW(p,±\t\)=A m 



N {1) - ifirfsn Nfa + Rfa - M + /i 5 r 3 R\ x) - Nfa - 2/i 



-z^ + Rf 1} - 2/v m5 r 2 - 2AT (1) ^ 5 r 3 ) (lC ± ^sinh£ (1) j + ((1) (2)), (17) 



with: 
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{1> '^ 4((M , h + l/af - /i 2 ) sinh£ (1) , (2) (cosh£ (1) , (2) - cosh£ (2))(1) ) ' ( ' 

JV( 1 ) j ( 2 ) = Mo, h + ~ (l - cosh S(i) j( 2) Y ^(i),( 2 ) = IC 2 + fi a + fi S - ^ sinh 2 £(i),( 2 ) , (19) 

1 3 2 3 

£ = - ^7jsin(pja), M 0>a; = m 0lX + - ^ sin 2 (p i a/2), x G {/, a}, (20) 

a i=i a i=i 

and the poles of the propagator in the energy-momentum space representation: 



, _ -b^pyb 2 -Aac , , 

cosh£ {1)i(2) = + V 2 . , (21) 

where: 

5 = ;KK + y 2 -4 (22) 

6 = -/i 2 (M ^ + -) --s(m ^ + -), (23) 

ay ay a \ ay 

c = s 2 - 4 / i 2 /i 2 - 4/i 2 (M 0th + ^ , (24) 

s=(M 0)h +-) + ^ + /C 2 + /i 2 + /£ (25) 

This analytical result for the heavy propagator has been checked by comparing with 
numerically computed propagators for various spatial lattice extensions and quark masses. 

The corresponding light twisted mass propagator has been calculated in [13]. It reads: 

S w (p,±|t|) = B^(l -cosh E + aMoj) ±^smkE - iaK - ia^r^, (26) 
where K. and M Q j are defined in (20) and 

_ e-^IA* _ (oMp,, + l) 2 + /C 2 a 2 + /i 2 a 2 + 1 

2sinh£(l + aM 0i )' 2(l + a¥ oi ) ' { ' 



4 Correlation matrices for K and D mesons 



To study cut-off effects at tree-level of perturbation theory for Nf = 2+1+1 quark flavours, 
we consider the spectrum of K and D mesons. As already discussed in Section 2.1.2, such 
mesons consist of the light up/down antiquark and either the heavy strange or the heavy 
charm quark. Since there is no gluonic field at tree-level, both quarks are free particles. 
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To create such mesons with well defined quantum numbers, we apply meson creation 
operators C(h,r) to the vacuum state \fl). In the physical basis in continuum QCD, a 
possible choice of appropriate operators is given by: 



rf 3 x^ (u) (x,t)rV' w (x,t). 



(28) 



The heavy flavour index h £ {s,c} determines whether the K meson or the D meson is 
created, the 4x4 matrix T G {75, 1} realizes the total angular momentum J = and 
either the negative parity V = — or the positive parity V = + and the integration over 
space / d 3 x yields the total momentum P = 0, i.e. assures that the light and the heavy 
quark have opposite momenta ±p. In principle, one could also fix the individual quark 
momenta ±p by including a second integration over space, but we prefer to consider not 
only the ground state, but also higher states, to have a situation, which is more like one 
of the interacting case, i.e. beyond tree-level, where one cannot get rid of excited states 
at the stage of operator construction, see Refs. [6,7] for an investigation of meson mass 
determinations in the interacting case. 

On the lattice using the twisted mass formalism an equivalent set of meson creation 
operators is given by: 



Note, however, that at finite lattice spacing, even after a rotation to the physical basis, 
parity and heavy flavour are only approximate quantum numbers, because of explicit 
twisted mass flavour and parity breaking. For a detailed discussion of these issues we refer 
to [6,7]. To extract meson masses, we first calculate 4x4 correlation matrices with the 
four operators (29): 



Details of this calculation, which uses the propagators from the previous section, are 
presented in App. B. Then, we extract the masses of the K and D mesons, as explained 
in details in App. C and [6]. 



5 Numerical results 

In this section, we will present some results for the continuum limit scaling of meson masses 
using maximally twisted mass fermions. In the investigation below, we will employ both 
the optimal and non-optimal definitions of maximal twist. 




(29) 



X 



C(/n,ri) I (/» a ,ra)(<i -^2) = (Q\(0( hltTl ){ti)^ {h2) r 2 )(t 2 )\ty 



(30) 
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5. 1 Setup and physical parameters 



We consider a setup, which is reminiscent of QCD with Nf = 2 + 1 + 1 flavours of quarks: 
there is the degenerate doublet of light fermions (masses m u ^) and the non-degenerate 
doublet of significantly heavier fermions (masses m s 22m n ^ and m c ~ 284m Uj d), where 
the ratios of masses have been chosen as observed in the nature for the up/down, strange 
and charm quarks. However, in contrast to QCD, we consider the fermions (to which we 
will also refer as "quarks") at tree-level, i.e. there are no interactions of any kind. 

We consider a 3-dimensional spatial volume of extension iV 3 with periodic boundary 
conditions. Having fixed the ratios of quark masses, there are two dimensionful parameters 
remaining, N and m u ^. Their dimensionless product Nm Ut d fully determines the physical 
situation. We choose Nm u ^ = 0.01. Roughly speaking, this assures that even for our 
smallest lattices, the heavy charm quark mass m c in lattice units is smaller than 1 (see 
below for more details). 

5.2 Continuum limit and lattice parameters 

For our lattice computations, the following set of parameters has to be chosen: the number 
of lattice sites N in spatial direction (i.e. the spatial lattice volume is N 3 ), the twisted 
quark masses fi q , fi a , fi$ (in units of N = 1/a) and the untwisted quark masses m 0j z and 
mo,/, (in units of N — 1/a). 

At tree-level of perturbation theory, the continuum limit (a — > 0) and the infinite volume 
limit (N — > oo) are equivalent. This means that the role of lattice spacing a is played 
simply by the inverse of the number of lattice points N. 

To recover in the continuum limit the setup described in the previous subsection, the 
lattice parameters have to be chosen in the following way, if the optimal definition of 
maximal twist is used: 

• Untwisted quark masses: 

m ,i = m 0th = 0. (31) 

• Twisted quark masses: 

Nfi q = m M , N(fi a - [is) = m s , N(fi a + /i s ) = m c . (32) 

When approaching the continuum limit, we set the twisted quark masses as specified in 
(32), i.e. we use this choice of twisted quark masses also at finite lattice spacing. The 
actual values of the strange and the charm quark mass for our numerical simulations are 
chosen such that they correspond to the ratios of central values of quark mass estimates 
published by the Particle Data Group [20] - m s /m u ^ = 21.5901 and m c /m ux i = 284.0910- 

1 Of course, since we are studying here an unphysical situation which, however, can illustrate 
the size of cut-off effects, the definite values for the quark masses used are not very important 
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This implies that since N(n a + fi$) = Nm c m 0.01 x 284 = 2.84, even for our smallest 
lattices, corresponding to N = 4, the heaviest quark mass in lattice units (equal to fia+Hs) 
is smaller than 1. 

Regarding the untwisted quark masses, we study different ways of approaching the contin- 
uum limit. A particular and optimal choice is m 0i / = rriQ ^ = 0. However, any other choice 
Nrrioj = 0(1/ N) and Nm Q)h = 0(1/ N) also corresponds to maximal twist. Therefore, we 
also study the choices: 



Nmo t i = Nm 0th = ^, (33) 

with q and q, arbitrary, but constant. Note, however, that in practice such choices should 
not be considered since -while still keeping the (9(a)-improvement of the theory- they 
might lead to the chirally enhanced cut-off effects, as discussed in Refs. [15,16,17,18,19]. 



5.3 Numerical study of the continuum limit 



We study the continuum limit by performing computations with various lattice volumes, 
ranging from N 3 = 4 3 to iV 3 = 96 3 . In Figs. 1, 2 and 4 we investigate the contin- 
uum limit of rriK and ttid for different realizations of maximal twist, c = q — Ch £ 
{0, 0.04, 0.06, 0.08, 0.1} (Fig. 1), c = q = c h G {0, 1.0, 1.5, 2.5} (Fig. 2) and q = 0, c h = 0.1 
or q = 0.1, Ch = or q = q, = 0.1 (Fig. 4). The curves in the plots correspond to fits of 
quartic polynomials in 1/N 2 : 

Nm n ^ D = a + a 1 ^ + a 2 ^^j + a 3 Q^) + a 4 (^) • (34) 

In the former case, i.e. for relatively small values of the parameter c and large lattices 
(N = 32 to N = 96), we observe a linear dependence in 1/N 2 = a 2 (i.e. the values of 
higher order coefficients a n (n > 2) in (34) are very small). This is, of course, expected, 
since Wilson twisted mass lattice fermions at maximal twist guarantee the absence of 
0(a) discretization effects. 

On the other hand, for large values of c = q = (as shown in Fig. 2), both K and D meson 
masses as functions of 1/N 2 clearly exhibit a non- vanishing curvature increasing with the 
value of c. However, this curvature can be well described by higher order corrections 
in 1/N 2 , as required from an (9(a)-improved theory. Nevertheless, when choosing the 
coefficients q^ large, the theory is more and more tuned towards the Wilson quark action 
showing large cut-off effects. For very large c values, also the extracted continuum limit 
values of meson masses are not reliable - because of the large curvature, in order to 
extrapolate to the continuum limit, one would need lattices with N > 96 in such case. Of 

and any reasonable value of the ratio of quark masses would be sufficient. Nevertheless, we have 
used the accurate values of ref. [20] to reflect the physical ratios of the quark masses as close as 
possible. 
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Fig. 1. The cut-off effects and the continuum iimit of: (a) kaon mass, (b) D meson mass, for 
different realizations of maximal twist c = q = Ch G {0,0.04,0.06,0.08,0.1} and 32 < N < 96 
lattices. 



N-m q =0.01, L-M K con =0.22591 , c=c,=c h N-m q =0.01, L-M D con '=2.85091 , c=c,=c h 




Fig. 2. The cut-off effects and the continuum limit of: (a) kaon mass, (b) D meson mass, for 
different realizations of maximal twist c = q = Ch G {0, 1.0, 1.5, 2.5} and 32 < N < 96 lattices. 

course, the values of q ^ used here are exceptionally large, leading to cases of non-optimal 
tuning to maximal twist that would never be used in practice. Hence, our investigation of 
these non-optimal tuning conditions serves solely illustrative purposes to indicate when 
higher order corrections in 1/N 2 become relevant. 

The values obtained for the ao coefficient correspond to meson masses extrapolated to the 
continuum limit. Within our numerical precision, they agree with the expected continuum 
results, i.e. the respective sums of the corresponding quark masses Nm u ^ = 0.010, Nm s = 
0.21591 and Nm c = 2.84091. 

The slope parameter a± describes the magnitude of 0(a 2 ) discretization effects. The ex- 
tracted values of this parameter are shown in Fig. 3 for the cases c = q = Ch and 
q = 0, Ch > 0. In the case of q = Ch = (Fig. 3(a)), the cut-off effects are significantly 
larger for the D meson than for the K meson, while the smallest ones are observed for the 
pion. Considering the relative discretization effects, given by the ratios ai/ao, we observe 
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Fig. 3. The slope a\ in eq. (34) for the pion, kaon and D meson masses vs. the values of q and 
Ch. (a) The case c = q = c/j. (b) The case q = 0, > 0. 
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Fig. 4. The cut-off effects and the continuum limit of: (a) kaon mass, (b) D meson mass, for 
different realizations of maximal twist: (q = 0, / 0), (q ^ 0, Ch = 0), (q = c/j / 0). 

that the cut-off effects are ~ 10 2 times larger for the D meson than for the K meson and 
~ 10 3 times larger for the kaon than for the pion. Moreover, the size of discretization 
effects increases with increasing values of the parameter c = q = Ch, but the sensitivity 
to this parameter is clearly the largest for the pion mass and it is very weak in the case 
of the D meson mass. This is an expected feature of the considered setup, since the non- 
vanishing bare Wilson quark mass, introduced by non-zero values of the parameter c, is 
relatively large in comparison with the pion mass and almost negligible as compared to 
the D meson mass (unless very large values of c are considered). The case q — 0, > 
(Fig. 3(b)) will be referred to below. 

However, as we have already mentioned, 0{a 2 m 2 q ) effects are expected to be present in 
0(a)-improved theories and they can become important in theories with heavy strange and 
charm quarks. Therefore, the 0(a 2 ) discretization effects given by the extracted parameter 
di contain the "pure" 0(a 2 ) effects and in addition the 0(a 2 m 2 ) effects. 

In order to disentangle 0(a 2 ) and 0(a 2 m 2 ) effects, we have computed the dependence of 
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Fig. 5. The 0(a 2 ), 0(a 2 m 2 ) and 0{a 2 m 4 ) effects (where q = s,c) in: (a) kaon mass, (b) D meson 
mass. The case of optimal tuning to maximal twist (q = Ch = 0). 
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Fig. 6. The extracted values of the parameters (a) &n, (b) b±, for different values of the parameter 
c = q = c/j. &o and b\ are defined by eq. (35). 

the K and D meson masses on Nm s and Nm c , respectively, for fixed lattice size N. Then, 
we have fitted the following function to the lattice data: 



NM lat - NM, 



cont 



b 1 + b {Nmq) 



(Nm q ) 4 
N 2 ' 



(35) 



where: NM\ & t - the lattice value of the K or D meson mass, NM cont = Nm Ut d + Nm q - 
its continuum (N — > oo) counterpart, q = sore. The dependences of the K and D meson 
masses on Nm q , together with the fit of the above functional form, are depicted in Fig. 5 
for the case q = q, = 0. The extracted values of bo, bi and b<i are shown in Table 1 for 
q = Ch = 0, q = Ch = 0.1 and q = 0, Ch = 0.1. For our fit, we have used TV = 48. 
However, we have checked for other lattice sizes (N = 16, 32, 64) that the values of bo, by 
and 62 remain the same, within numerical precision. 

Our fits are summarized in Fig. 6 (for the case c = q = qj and Fig. 7 (for q = 0, > 0). 
Fig. 6(a) shows the magnitude of "pure" 0(a 2 ) effects for the pion, the kaon and the D 
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N-m q =0.01, 




N-m q =0.01, C|=0 



0.3 0.4 0.5 0.6 0.7 




(a) 



(b) 



Fig. 7. The extracted values of the parameters (a) bo, (b) b±, for different values of the parameter 
Ch, with q = 0. &o and b\ are defined by eq. (35). 





c = q = 


~-c h = 


c = q = 


c h = 0.1 


q = 0, c h = 0.1 


Coefficient 


K meson 


D meson 


K meson 


D meson 


K meson 


D meson 




0.000571 


-0.552 


0.499 


-0.199 


0.00442 


-0.596 


h 


-0.0373 


0.214 


-0.0848 


0.198 


-0.0848 


0.183 


b 2 


-0.159 


0.00540 


-0.0623 


0.00550 


-0.0623 


0.00592 



Table 1 

The fitting coefficients of eq. (35). The fits are shown in Fig. 5 for the case N 



48, Q = c h = 0. 





c = q = c h = 


c = ci = c h = 0.1 


q = 0, c h = 0.1 


Effects 


K meson 


D meson 


K meson 


D meson 


K meson 


D meson 


0(o a ) 
0(a 2 m 2 ) 
0{a 2 m]) 


2.5 • 10- y 
-7.6 • 10~ 7 
-1.5 • 10~ 7 


-2.4- 10" 4 
7.5 • 10" 4 
1.5 • 10" 4 


2164 • 10" 7 
-17- 10~ 7 
-0.6 • 10~ 7 


-0.9 • 10" 4 
6.9 • 10" 4 
1.6 • 10" 4 


19.2 • 10~ 7 
-17.2 • 10~ 7 
-0.6 • 10~ 7 


-2.6 • 10" 4 
6.4 • 10" 4 
1.7 • 10~ 4 


sum 


-6.6 • 10" Y 


6.6 • 10~ 4 


2146 • 10- y 


7.6 • 10" 4 


1.4 ■ 10" 7 


5.5 • 10~ 4 



Table 2 

The decomposition of the difference A^M[ at — iVM cont for the K and D mesons. The lattice size 
N = 48, the strange quark mass Nm s = 0.21591 and the charm quark mass Nm c = 2.84091. 

meson. In the region of small values of the parameter c, these effects are the largest for 
the D meson and comparable to each other for the pion and the kaon. Moreover, the curve 
b (c) coincides with the curve ai(c) in Fig. 3(a), since in the case of the pion C(a 2 m^ d ) 
effects are negligible, due to the small value of light quark masses. We also observe that 
the value of bo increases with increasing c. However, since bo < for the D meson in this 
region, the size of 0(a 2 ) effects decreases when the parameter c increases, which means 
that the effects of non-optimal tuning can partially cancel 0{a 2 ) effects. For a particular 
value of the parameter c (around 0.13) a complete cancellation can even occur (i.e b = 
can result). 

Fig. 6(b) shows the magnitude of 0(a 2 m 2 ) effects for the kaon and the D meson. These 
effects depend only slightly on the value of the parameter c. Interestingly, the sign of 
0(a 2 m 2 ) is opposite to the sign of the "pure" 0(a 2 ) effects, which means that a partial 
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cancellation between these two types of effects occurs. 

Further insight into the role of a method of tuning to maximal twist can be obtained by 
analyzing the dependence of cut-off effects for q = and Ch > 0, i.e. if the Wilson mass 
mo,/ is precisely tuned to maximal twist, while the Wilson mass mo } h is non-optimally 
tuned (Fig. 7). In such case, the 0(a 2 ) effects are again much larger for the D meson than 
for the kaon. What is more, the magnitude of these effects is much smaller for the kaon 
than in the case q = c^. This implies that non-optimal tuning in the light sector leads to 
much larger effects than non-optimal tuning in the heavy sector. Still, both effects tend to 
increase the value of the coefficient bo- The situation is different for the D meson. There, 
the effects of non-optimal tuning of m ,/ and mo,/, have opposite signs and thus the "pure" 
0(a 2 ) effects in the case q = 0, q^ > increase with increasing value of the parameter c^. 
In this way, the partial or complete cancellation of 0(a 2 ) effects that we have observed 
in the case q = Ch can be attributed to non-optimal tuning of mo,; and not mo,h. 

Regarding the size of 0(a 2 m 2 ) effects in the case q = 0, Ch > 0, Fig. 7(b) shows that 
again non-optimal tuning in the light and in the heavy sector can have opposite effects. In 
the case of the D meson, a partial cancellation of 0(a 2 m 2 ) effects can occur for positive 
values of the parameter Ch, but in the case of the kaon, the effects of non-optimal tuning 
of the Wilson mass m 0i/l reinforce the magnitude of 0(a 2 m 2 ) effects. 

It is also interesting to see the relative contribution of 0(a 2 ), 0(a 2 m 2 ) and 0{a 2 m^) 
discretization effects to the difference NM\ at — NM cont for some chosen values of the 
lattice size and the strange and charm quark masses. Here we choose again N = 48 and the 
quark masses considered in the earlier tests described in this section, i.e. Nm s = 0.21591, 
Nm c = 2.84091. The decomposition of the difference NMi at — NM cont for such parameter 
values is shown in Table 2, for the cases q = = 0, q = Ch = 0.1 and q = 0, = 0.1. 

The size of 0(a 2 ), 0(a 2 m 2 ) and 0(a 2 m q ) cut-off effects is of the same order of magnitude 
in the case of optimal tuning to maximal twist in both the light and the heavy sector. 
However, in this case the largest contribution to the difference NM\^ — NM cont is the one 
of 0(a 2 m 2 ) effects, for both the K and D meson. Moreover, the relative sign of 0(a 2 ) and 
0(a 2 m 2 ) effects is different. The 0{a 2 rn^) effects are roughly a factor of 5 smaller than 
0(a 2 m 2 ) effects in the case of both meson masses. Again, the overall size of discretization 
effects is much more important for the D meson than for the kaon, both if we consider 
absolute and relative cut-off effects. 

In the case of non-optimal tuning of both mo,/ and mo,/i, the 0(a 2 ) effects (induced by non- 
optimal tuning particularly in the light sector) are much more important than 0(a 2 m 2 ) 
effects in the kaon case, but the latter still dominate in the case of the D meson. 

An interesting situation occurs in the case q = 0, Ch = 0.1 for the kaon, where we observe 
that 0(a 2 ) effects are almost equal in magnitude, but of opposite sign to the 0(a 2 m 2 ) 
effects. Thus, an almost exact cancellation occurs and the overall size of cut-off effects 
is much smaller than in the case of optimal tuning in both sectors. Such decrease of the 
overall size of discretization effects with respect to the optimal tuning case occurs also for 
the D meson, which is the effect that can be clearly observed in Fig. 3(b). 

To summarize, we observe a very intricate interplay of different types of effects - 0(a 2 ), 
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0(a 2 m 2 ) and even 0(a 2 rn^) effects can become sizable. A partial or complete cancellation 
between 0(a 2 ) and 0(a 2 m 2 ) effects is possible and moreover such cancellation can also 
occur between effects of non-optimal tuning to maximal twist in the light and heavy 
sectors. Needless to say, the behaviour in the interacting theory should be expected to be 
even more complex. 



6 Conclusions 



In this paper, we have provided an analytical basis for studying lattice spacing effects at 
tree-level of perturbation theory for maximally twisted mass Wilson quarks, when both 
the heavy quark doublet and the light one are included. Particularly, we have calculated 
the 4-dimensional momentum space and time-momentum frame quark propagators in the 
heavy sector and constructed the matrix of correlation functions for the K and D mesons. 

We have investigated the scaling of the kaon and the D meson masses with the lattice 
spacing, for optimal and non-optimal tuning to maximal twist. We have clearly verified 
that the lattice spacing effects appear in even powers of 1/N = a, as expected from 
the general automatic (9(a)-improvement of maximally twisted mass lattice QCD. More 
precisely, for the case of optimal tuning (given by the introduced parameter c = q = Ch = 
0, where q and Ch determine the light and heavy Wilson mass: Nm j = q/N, iVm 0) h = 
Ch/N) and non-optimal tuning with small c values, we have observed the linear dependence 
of meson masses in 1/N 2 = a 2 , while for sufficiently large c values, the 0{a i ) and higher 
order corrections can become important. 

We have also disentangled the 0(a 2 ), 0(a 2 m 2 ) and 0{a 2 rn^) discretization effects (where 
m q is the mass of the strange or the charm quark) in the kaon and the D meson masses and 
we have found a multitude of competing effects. The overall size of discretization effects 
can considerably depend on the details of tuning to maximal twist in both the light and the 
heavy sector. Partial cancellations (or even full cancellations for some particular values 
of parameters) can occur between the 0(a 2 ) effects and the 0(a 2 m 2 ) effects and even 
between effects of non-optimal tuning in the light sector and the effects of non-optimal 
tuning in the heavy sector. Thus, non-optimal tuning can in certain instances decrease 
the overall size of discretization effects. 

We believe that although our results do not provide any proofs for the interacting case, 
they give a clear warning: one should expect a very complex interplay of different effects 
and thus when physical observables in the strange and especially charm sectors are evalu- 
ated in full QCD, the enhanced cut-off effects should carefully be taken into account and 
physical results obtained at a single value of the lattice spacing can differ significantly 
from the continuum results. 
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Appendix A 



The kernel of the heavy quark propagator: 



(36) 



n=i 



+ 



x = na, y = ma, n, m e Z. 

Further, we skip colour, flavour and Lorentz indices. 

After discrete Fourier transformation: 



OxyOABi 



K{p) = Y J K{x-y)e 



-ip(x-y) 



(37) 



x-y 



we get: 



k (p) = - X) 7i s m (ft«) + - (1 - cos(p 4 a)) + sin2 Hr ) + m °> h + ^/V/sTl + Hst 3 - (38 



a i=i 



Denote: 

Af ,h = mo,h + - > sm 

a ti V 2 



2 / Pi a 



N = M 0>h {p) + -{l-cos{p 4 a)), 

(X 



1 3 1 

JC = - ^2 ji sin(pia), R 2 = fij + fil + K, 2 + — sm 2 (p 4 a) 



i=i 



Then, the matrix operator K(p) takes the form: 

= N + fl S T 3 + + i/I a l5n- 

The propagator can be calculated from the following relation: 

K\p) 



S^ h \p)K{p) = 1, S™(p) 



K(p)&(py 



(39) 
(40) 

(41) 
(42) 
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Using the properties of gamma matrices, we get the heavy twisted mass propagator: 

n(h) ( , = N + ^3 ~ * gU ^ sin^q) - tfi^n 

W N2 + R* + 2Nvst 3 + 2^[j, 615 t 2 ■ 1 } 

The denominator has the form n + mT 3 + lr 2 , where n, m and I are numbers. To eliminate 
the flavour structure from the denominator, we rearrange: 

(n + mr 3 + lr 2 )(n - mr 3 - lr 2 ) — (n + mr 3 )(n - mr 3 ) + W 2 (n - mr 3 ) - (n + mT 3 )lr 2 -I 2 = 



n 2 + mnr 3 — nmr 3 — m 2 + lnr 2 — lmr 2 T 3 — nlr 2 — mlr 3 T 2 — I 2 = n 2 — m 2 — I 2 . (44) 
After this, we obtain: 

{h)( v. _ (N + fi S T 3 - i g^=i y sin(p M q) - i/yysTi) (N 2 + R 2 - 2Nfi 5 r 3 - 2/i CT /i 3 7 5 r 2 ) 
^ ~ (iV 2 + P 2 ) 2 - 4iV 2 /^ - 474a*I ' 

(45) 

The denominator of the last expression has two zeros which correspond to the poles of 
the propagator: 



-61 =F v /(W) 2 -4(al)(cl) 
coshP li2 = cos(p 4 a) li2 = . (46) 

We calculate the residues of the function: 

/ = 7 2 x ^ (47) 

' (M ,h + i) - /i 2 J (cosh E - cosh P x ) (cosh £ - cosh E 2 ) 

The first residue is taken at the point E = E\. 



, -s ia p t Pi (Pi) , Ar 

(res f) E=El = — e"^ — , (4* 

^o,/i + 1 ) — A*| ) sinh Pi (cosh E 1 - cosh E 2 ) 



where: 

Pi(P0 = N^N 2 + R\- 2/i 2 ) + N (Rj - N 2 - 2/i 2 )r 3 - i(N 2 + R\) (jC + *3± sinh P x ) + 



22/1^/15 ^/C + ^ sinh P^ 75 r 2 + 2iN x ^ 5 \ K. + ^ sinh P^ r 3 - i[i a (N 2 + R\- 2/i 2 )7 5 Ti, 

(49) 

iVi = M 0A + -(1 -cosh Ei), Rj ^ ]C 2 + ^ + ^ - ^ sinh 2 El (50) 
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The second residue is at the point E = E 2 : 



(res f) E=E2 = — e fia »y- - 2 r , (51) 

' [M ,h + J) - A*i ) s inh ^2 (cosh £ 2 - cosh Si) 



where: 

Pi^a) = iV 2 (iV 2 2 + i? 2 - 2^ 2 ) + ^(i? 2 - N% - 2/4)r 3 - f(iV| + i^) f/C + ^ si 



2ifi a fi S (jC + sinh S 2 ^ 75 r 2 + 2iN 2 fi s (fC + sinh S 2 ^ r 3 - 2/i CT (iV 2 2 + Rj - 2fi 2 s )^ 5 r u 

(52) 

^2 = M , h + -(1 -cosh£ 2 ), R\ = K 2 + \i a + /i 5 - -j sinh 2 E 2 . (53) 
a er 

The propagator is then the sum: 

Soc(p, t) = 2%i [(res f) E =E 1 + (res f) E =E 2 ] ■ (54) 
Appendix B 

The elements of the correlation function matrix are calculated for infinite time using 
the expressions for infinite time twisted mass propagators (17) and (26). We show the 
calculation of one of them. We use the Fourier transformation, the definition of the 5- 
function and gamma algebra: 



Cn(t) = E(^ (u) W75^ (s) W) t (^ (n) (0)75^ {s) (0)) 



E Tr{ 75 S fm (x, t; 0, 0) 75 S Ss (0, 0; x, t)} 



,L 3 

x \ p 



(55) 



1 E Tr f J3 E E e^-^S™(p, t) l5 S Ss (q, -t) 75 



L 3 



5 x 



^ E Tr ^ E <7)^(£ *) 75 ^(g, -t) 7 5 
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= J? E Tr P**® *)75^ S (P, -t)7s) • 
p 

The components of the light and heavy quark propagators can be written in the following 
way: 

3 3 

S~ SS (P, -t) = E 7/, + 74/4 + 1 • h, S nu (p, *) = £ lidi + 74^ + 7s4 + 1-4- (56) 

i=l t=l 

Thus, we get: 



Tr {S^(p, t) 75 S 5s (p, -t) lB } = Tr{ £ "4/, - 4/4 + 4/e). 

^ i=i J 



(57) 



We identify the corresponding coefficients fi and di, where % = 1, .., 6 in the expressions 
(56) and (57) and obtain: 

C u (t) = £ {aCXV^ 2 + #D + 5(iV 2 2 + Rl) + 2/x 5 (AiV 1 + £?iV 2 )] + (58) 

v 

Csinh^ 
+ - 



A(N{ + Rf) sinh E 1 + £(iV 2 + i^) sinh E 2 + 2pi 5 (AN 1 sinh E 1 + 5iV 2 sinh £ 2 ) + 
+C(1 - coshEj + aAfo.OfAWi^ + i? 2 - 2/$ + 5iV 2 (iV 2 + R\ - 2$)- 
-lxs[A(Rl - iV 2 - 2^ 2 ) + B{R\ - Nl - 2/4)]]). 



C 22 (t) = E((^ (U) (^)75^ (C) (^))¥ W (0)75^ (C) (0)) 



N c N d 



£ {aC7C 2 [A(iV 2 + i? 2 ) + 5(iV 2 + R 2 2 ) - 2fi s (AN! + BN 2 )]+ (59) 
p 

Csinh 



A(Nf + Rf) sinh £ x + B(N 2 2 + i? 2 ) sinh E 2 - 2fi S (AN 1 sinh ^ + BN 2 sinh £ 2 ) + 

a L 

+C(1 - coshEj + aM^AN^Nl + R\- 2/4) + BN 2 {Nj + R 2 2 - 2/4)+ 
+fi S [A(Rl - Ni - 2/i 2 ) + B{R\ - N\ - 2/, 2 )]]}, 

C33W = E((^ (u) (^)^ (s) (^))¥ (?l) (o)^ (s) (o)) = 
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NcNd 
L 3 

v 



{aCK?[A{Nl + R\) + B(N% + R 2 2 ) + 2 l i s {AN 1 + BN 2 )}+ (60) 

p 

Csinh E[ 



+ A(N{ + sinh E x + £(7V 2 + R z 2 ) sinh £ 2 + 2// 5 (ANi sinh E x + 5iV 2 sinh E 2 ) 

-C(l - coshEj + aM Qt i)[AN x (Nl + R\- 2/4) + BN 2 {Nj + R\ - 2/4)- 
-»s[A(Rl - Nl - 2/4) + B(Rj - Nl - 2/i 2 )]]}, 

c M (t) = x:((^w(x)^ (c) (x))^ (u) (o)^ (c) (o)) = 

.? 

= E {oC£ 2 [4(iV? + # 2 ) + fi(iV 2 2 + i? 2 2 ) - 2 yU5 (AiV 1 + BN 2 )]+ (61) 

p 

+ CsmhE 1 2 + - ^ ^ + B 2 + „ ginh ^ _ 2^ M JVi sinh ^ + £JV 2 sinh £ 2 ) 
a L 

-C(l - coshEj + aM^AN^Nl + R\- 2$) + BN 2 (N 2 + R 2 2 - 2/4)+ 
+pi s [A(Rl - Nl - 2/4) + B(Rl - Nl - 2/4)]]}, 



C48(t) = E((^ (lt) (^)^ w (^)) t ^ H (o)^ w (o)) = 

X 

^ aCfi qf i a {A(N* + R\- 2/i 2 ) + I?^ 2 + R\ - 2/4)}, (62) 



N c N d 



L 3 

p 



Cu(t) = E((^ W (^)^ (c) (^))¥ (u) (0)c(0)) 

X 

N c N d 



L 3 

p 



£ aC^M^iV 2 + ^ " 2/i 2 ) + B(N 2 + R 2 2 - 2/i 2 )}, (63) 



Csi(t) = E((^ (lt) (^)^ w (^)) t ^ H (o)75^ (s) (o)) 



Y l (iaCfji q ){AN 1 {Nf + R\- 2/i 2 ) + 5iV 2 (iV 2 + # 2 - 2/i 2 )- (64) 



_ N c Nd ^ 
p 

-» S [A(R 2 ~ Nl - 2/i 2 ) + B(R 2 2 - Nl - 2/i 2 )]}, 
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CisW = E((^ w (^75^ (s) (^)) t «(o)^ (s) (o)) 



= Y.(- iaC N){ANi{Nl + R\- 2/4) + BN 2 (Nl + R\ - 2/4)- 

p 

-H & [A{Rl - Nl - 2/4) + B(Rj - N 2 2 - 2/4)]}, 

= E((^ (lt) (^)^ (c) (^))¥ (u) (o)75^ (c) (o)) = 

.? 

2(i^ ff ){^JVi(JV? + i? 2 - 2/4) + £iV 2 (iV 2 + i? 2 - 2/4)+ 



N c N d 



L 3 

p 



+» s [A(Rl - iV 2 - 2/4) + B{R\ - Nl - 2/4)]}, 



= ^jr ^{-iaC N ){AN x {Nl + R\- 2/4) + BN 2 (Nl + R\ - 2/4) 
p 

+/i 5 [A( J R 2 - N\ - 2/4) + B(Rj - N 2 2 - 2/4)]}, 



= E((^ W (^)75^ (S) (^))¥ W (0)75^ (C) (0)) 

..r 



L 3 



£(-aC7/tyO {^(iV 2 + i? 2 - 2/4) + 5(iV 2 + i? 2 - 2/4)} 



CaiW = E((^ W (^)75^ (c) (x))^W(0)75^ (s) (0)) = 

X 

Vd E(-aC/i g /i CT ) {A(iV 2 + i? 2 - 2/4) + 5(iV 2 + R\ - 2/4)} , 



L 3 

v 



c u (t) = E((^ (lt) (^)75^ (s) (^))¥ (u) (o)V' (c) (o)) 



N c N d ^, _ „ Q , . „, 2zC 



E{-2«Ca/C 2 /i (7 // ( 5(A + 5) sinh E[^ S {A sinh E t + B sinh £ 2 ) 



L 3 ^ L ^ ov ' y a 

p 
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-iC/v(l - cosh£{ + aM^)[A{Nl + R\ — 2/x|) + 5(iVf + R\ - 2/4)}}, 



c 4 i(t) = E((^ (lt) (^)^ (c) (^))¥ (u) (o)75^ w (o)) 



NcNd^ „o , . 2iC 

L 3 - 



^2{2iCaK. 2 /j, tT fj. s (A + B) + sinh E[^ S (A sinh ^ + 5 sinh £ 2 ) + (71) 



+z'C/i CT (l - cosh£{ + aM ,i)[^(iV 1 2 + -Rf — 2/i 2 ) + 5(iV| + R\ - 2^)]}, 



c^W = E((^ (u) (^)^ (s) (^))¥ w (o)75V' (c) (o)) = 

= Y {-2iCaJC 2 fi a ^s(A + B)-— sinh E[fi a fi s (A sinh ^ + S sinh £ 2 ) + (72) 

p 

+iC/v(l - coshEj + aM ,)[A(iV 1 2 + R\ - 2/4) + 5(iV| + R\ - 2/x 2 )]} 



C 23 (t) = E((^ W (^)75^ (c) (^)) t ^ W (0)^ W (0)} = 

X 

= ^{2iCa/C 2 /i CT /i 5 (A + 5) + — sinh E\\i a \i&{A sinh £i + B sinh £ 2 )- (73) 

-iC/x CT (l - cosh£{ + aM o ,0[A(A r i + R\ - 2/i 2 ) + B(N% + R\ - 2/i 2 -)]}. 



Appendix C 

From the 4x4 correlation matrix in the twisted basis presented in App.B we get the 
masses of K and D mesons by solving a generalized eigenvalue problem [21]: 

]TC^(t)4 n) (Mo) = Y / C 3k (t )v i :\t,t )X^(t,t ), (74) 

k k 

where k runs over (h, T), h G {s, c} and T G {75, 1}, n = 0, 3. 

Eq. (74) leads to the following equation, giving the four effective masses v%g for the 
lattice with temporal extension T: 

\{n)U t() \ e -m ( ^\t,t )t _|_ e -m$(t,to)(T-t) 

(75) 



\( n )(t + 1, t ) e -^iff(*.<o )(t+l) + e -m^(t,t )(T-(t+l)) ' 
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Since we have performed the integration over p4, our masses correspond to the limit 
of infinite T. We consider finite t values, so the formula (75) becomes a ratio of two 
exponentials. For our tree level calculations, we have used t ^ 10. We have explored the 
dependence of the K and D masses on time and we have found that for the considered 
spatial extensions L the mass plateaus are reached at t < 100. 

We extract the meson masses from the plateaus at t G [50, 120] (depending on the spatial 
extension L of the lattices). In this range of t, our numerical precision gives plateaus of a 
very good quality and thus reliable data for the masses. 

To determine the parity and flavour content of the four effective masses we perform the 
approximate rotation to pseudo physical basis [6]. After the rotation, the correlation 
matrix takes the form: 

C$ r) (t) = M(u h u h )Cl r) (t)M(u h u h )l (76) 
The twist rotation matrix is orthogonal at maximal twist and has the form: 



M(uji, LO h ) 



( cos f cos 
- sin f sin f 
—i sin y cos f 7j L 

\-i cos f sin <f 



- sin f sin 

cos f cos <f 
—i cos y sin ^ 
-% sin f cos <f 



-% sin f cos 
-i cos f sin 
cos f cos f 
- sin f sin <f 



-i cos f sin 
-i sin f cos 
- sin T sin ^ 

cos f cos ^ 



(77) 



and connects the operator C^p) i R pseudo physical basis with the 0* h r ^ in the twisted 
mass one: 



<D = M{u u u h )Ol^ (Og r) )t = (^ r) )tM T (^,^) 
These bilinear operators are chosen in the following way: 



(7? 



(x {u hsx s \ 

X {u) lbX 

x (u) x s 
V x {u) x c J 



(79) 
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